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Abstract 

In order to study the relation between backbone and side chain ordering in proteins, we have 
performed multicanonical simulations of deka-peptide chains with various side groups. GIuiq, 
GlriiQ, AspiQ, AsriiQ, and Lysio were selected to cover a wide variety of possible interactions between 
the side chains of the monomers. All homopolymers undergo helix-coil transitions. We found that 
peptides with long side chains that are capable of hydrogen bonding, i.e. GIuiq, and GItliq, exhibit 
a second transition at lower temperatures connected with side chain ordering. This occurs in 
gas phase as well as in solvent, although the character of the side chain structure is different in 
each case. However, in polymers with short side chains capable of hydrogen bonding, i.e. Aspio 
and AsriiQ, side chain ordering takes place over a wide temperature range and exhibits no phase 
transition like character. Moreover, non-backbone hydrogen bonds show enhanced formation and 
fluctuations already at the helix-coil transition temperature, indicating competition between side 
chain and backbone hydrogen bond formation. Again, these results are qualitatively independent 
of the environment. Side chain ordering in LysiQ, whose side groups are long and polar, also 
takes place over a wide temperature range and exhibits no phase transition like character in both 
environments. Reasons for the observed chain length threshold and consequences from these results 
for protein folding are discussed. 
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I. INTRODUCTION 



In the last two decades, energy landscape and folding funnel paradign>i have led to an 
emerging understanding of the folding process in proteins. However, these concepts describe 
only the general characteristics of folding. Many details still remain poorly understood. 
One aspect is the role of side chain ordering in the folding process. Are side chain and 
backbone ordering coupled? Do the processes occur in a particular order? In principle, 
these questions can be investigated in silico. As one cannot employ lattice proteins and other 
minimalist models that ignore side chains, one has to rely on all-atom models of proteins. 
However, simulations of such protein models are often hampered by poor convergence^, and 
in particular slow side chain dynamics can pose problems^. Only with the development of 
generalized ensemble techniques^ such as parallel tempering^'^'^ or multicanonical sampling^"^ 
have all-atom simulations of small proteins (with up to ~ 50 residuesiS) become possible. 
In the present work we employed multicanonical sampling^ which was first introduced to 
protein science in Ref. lH 



In a previous publication^^ we reported results on backbone and side chain ordering 
in polyglutamic acid, GIuiq. Our results showed that - upon continuously lowering the 
temperature - GIuiq "folding" is a two-step process, starting with the secondary structure 
formation. Only after the backbone geometry is fixed the side chains order themselves at a 
much lower temperature. This scenario applied to the polymer in gas phase as well as in a 
solvent. It is quite remarkable that the scenario is independent of the particular environment 
since the side chain ordering had a different character in each case. In gas phase, side chains 
align themselves along the helical cylinder, stabilizing themselves by forming hydrogen bonds 
with each other, while in solvent the side chains extend into the solvent which screens them 
from forming hydrogen bonds. 

These findings for GIuiq, decoupling of backbone and side chain ordering independent 
of the environment, immediatly raise the question whether such a scenario is a common 
characteristic in protein folding, at least for hydrogen bond forming side chains. 

In order to test this hypothesis we are investigating in this contribution the dependence 
of side chain ordering on chain size and chemical properties. Glutamine (Gin), for example, 
is about the size of glutamic acid and is also able to participate in hydrogen bonds. Aspartic 
acid (Asp) and Asparagine (Asn) are also able to form hydrogen bonds, but have a smaller 
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size, while Lysine (Lys) has a larger polar side chain. Consequently, we have performed 
multicanonical simulations of Glnio, Aspio, ^-snio, and Lysio, and will compare them to 
extended results for GIuiq. 

Again, the folding scenarios found are independent of the environment, despite the fact 
that the side chain structure that arises is different for gas phase and in solvent. As an aside 
we note that the structural features found here confirm early investigations performed by 
H. Scheraga's group some 40 years agcM'^^'^^^. 

However, we observe a decoupling of backbone and side chain ordering transitions - as 
described above - only for GIuiq and GIuiq, i.e. for long hydrogen bond forming side chains. 
Polymers with short hydrogen bond forming side chains, i.e. AspiQ and Asuiq, exhibit no 
such decoupling. Instead, side chains already participate at the helix coil transition as 
competitors for the backbone hydrogen bonds. Furthermore, the final side chain ordering in 
these polymers takes place over a wide temperature range and exhibits no phase transition 
like character. Lysio whose side chain is long and polar follows the same szenario as Aspio 
and Asriio, albeit without the side chains competing at the helix-coil transition. 

II. METHODS 

Our simulations utilize the ECEPP/3 force field~ as implemented in the 2005 version 
of the program package SMMP— Here the interactions between the atoms within the 
homopolymer chain are approximated by a sum Eecepp/s consisting of electrostatic energy 
Ec, a Lennard- Jones term E^j, hydrogen-bonding term Ehb and a torsion energy Etot- 
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where is the distance between the atoms i and j, C,i is the l-th torsion angle, and energies 
are measured in Kcal/mol. The protein-solvent interactions are approximated by a solvent 



-Secepp/s — Ec + El J + Ehb + Exor 
_ ^ 332qiqj 





(1) 



3 



accessible surface term 

Esoiv = '^o-iAi . (2) 

i 

The sum goes over the solvent accessible areas Ai of all atoms i weighted by solvation 
parameters cTj as determined in Ref. |20|, a common choice when the ECEPP/3 force field 
is utilized. Our previous experiences^i^^ have shown that Egoiv reproduces the effects of 
protein- water interaction qualitatively correct. However, thetemperature scale is often dis- 
torted, leading, for instance, to transitions at temperatures where in nature water would be 
vaporized. This problem can be remedied, however, by renormalization of the temperature 
scale upon comparison with experiments. 

The above defined energy function leads to a landscape that is characterized by a mul- 
titude of minima separated by high barriers. As the probability to cross an energy barrier 
of height AE is given byexp(— AE'/Zc^T), ks being the Boltzmann constant, it follows that 
extremely long runs are necessary to obtain sufficient statistics in regular canonical simula- 
tions at low temperatures. Hence, in order to enhance sampling we rely on the multicanonical 
approach^'^ asdescribed in Ref. [llj. Here, configurations are weighted with a term wmu{E) 
determined iteratively such that the probability distribution obeys 

Pmu{E) oc n{E)wMu{E) ^ const , (3) 

where n{E) is the spectral density of the system. Thermodynamic averages of an observable 
< O > at temperature T are obtained by re-weighting^: 

JdxO{x)e-^(^y'-^/wMu[E{x)] 

where x counts the configurations of the system. 

After determining the multicanonical weights Wmu{E) we have performed multicanoni- 
calsimulations of 5 x 10® sweeps. Each sweep consists of Nf Metropolis steps that try to 
updateeach of the Nj dihedral angles (the degrees of freedom in our system) once. Here, 
Nf = 70,70,60,60,80 for Gluu), Glnio, Aspio, Astiiq, and Lysio, respectively. Every 10 
sweeps various quantities are measured and written to a file for further analysis. These 
includethe energy E with its respective contributions from Eq. ([T]) and - in the case of the 
simulations in solvent - from the protein-solvent interaction energy Egoiv The radius of 
gyration Vgy as a measure of the geometrical size, and the number of helical residues nfj, 
i.e. residues where the pair of dihedral angles {(pyip) takes values in the range (—70° ± 30°, 



—37° ±30°)^. Finally, we monitor the total number of hydrogen bonds, heb- Note that hy- 
drogen bonds along the backbone span four residues in the sequence. We therefore monitor, 
in addition, the number of hydrogen bonds between residues that are closer in the sequence, 
denoted by nfj^. This number is a lower limit to the number of non-backbone hydrogen 
bonds^i. 

III. RESULTS AND DISCUSSIONS 

Our aim is to study through multicanonical sampling the relationship between side chain 
ordering and other transitions for five specifically selected homopolymers. We are interested 
particularly in the behavior of these molecules at low temperatures where their energies 
are low, too. Hence, a measure for the quality of our results is the number of independent 
low-energy states. As a multicanonical simulation performs a random walk in energy space, 
see Fig. [H a lower limit for this number of independent states can be obtained from the 
number ut of passages from a suitably defined low-energy state to high-energies, and back, 
sometimes also called "tunneling processes". Since we wanted to obtain information about 
transitions that occur in the temperature range 2007^ to 600i^', we choose the average 
energies at T = I'bQK and at T = 750/^ as reference energies in each case. These reference 
energies and the numbers of tunneling processes for each simulation are reported in Table [B 
As in previous investigations, we found rix > 100 to be sufficent for a good quality of the 
results. As a check we performed a 10'' sweep simulation of Gluio in vacuum while all other 
simulations were 5 x 10^ sweeps. 

We first investigate the case of molecules in gas phase. Fig. [2] displays the specific heat 
per molecule. 



as a function of temperature for all five polymers. In each case one observes a peak at a 
temperature Ti in the range 450^ to 600K with quite large half- widths between 150i^ and 
200A'. The individual peak temperatures Ti and other properties are listed in Table [Tll 

The corresponding plot of the helicity in Fig. [3] shows that these peaks separate a high 
temperature region where the backbone has no ordering from a region where temperatures 
are low enough to allow the formation of an a-helix. As one can see from the monotonic drop 
in the average radius of gyration Vgy that is shown in the inlay, the helix-coil transition is 




(5) 
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TABLE I: Number of tunneling processes, n^, between the average energy at T = 150K, Ei^qx, 
and the average energy at T = 750K, Ei^qk in each simulation (5 x 10^ sweeps). The number in 
parenthesis is the error in the last digit. 





riT 


vacuum 
Ei^QK [Kcal/mol] 


^750X [Kcal/mol] 


riT 


solvent 
Ei^m [Kcal/mol] 


E750K [Kcal/mol] 




303 


-132.4(2) 


-42.6(1) 


183 


-210.08(5) 


-124.61(8) 


Glnio 


186 


-161.3(3) 


-76.29(8) 


251 


-232.42(5) 


-154.10(4) 


Lysio 


302 


-33.6(1) 


55.53(4) 


149 


-100.23(5) 


-10.24(6) 


AsTliQ 


168 


-195.65(10) 


-121.14(5) 


215 


-267.52(4) 


-196.86(4) 




118 


-169.7(1) 


-90.8(1) 


146 


-243.2(9) 


-169.6(9) 



"The GluiQ simulation in vacuum was 10^ sweeps for testing reasons. 



also connected with a collapse of the molecule. Below the transition Tgy stabilizes, reflecting 
the stable helical structure that has been reached. 

The transition temperatures are largest for Gluio and Glriio, and lowest for Asriio and 
AspiQ, the individual values being very close to each other in each block. The transition 
temperature for LysiQ is intermediate between both blocks. Since the side chains of Gluio 
and Glriio are larger than in the cases of Asriio and Aspio, they provide sterical hindrances 
to backbone conformations, leading to a decrease of the backbone entropy. As the transition 
is driven by entropy, and the energetic properties of these four heteropolmers are similar, 
this leads to a higher transition temperature for the polymers with larger side chains. This 
argument also explains the relationship of the transition temperature for LysiQ to that of 
a homopolymer of similar energetics but much shorter side chain, Alaio with a transition 



temperature Ti = 427i^', see Ref. 
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and Table [TTl 

However, this entropic argument does not suffice to understand why the transition tem- 
perature of Lysio is still smaller than that of GIuiq and GIuiq, although Lysio has even longer 
side chains than those two molecules. In order to explain this feature we have to take the 
different energetics into account, i.e. the fact that the side chains of GIuiq and GIuiq can 
participate in hydrogen bonds, which are stronger than it would be for the case of Lysio- 
These additional energy contributions shift the transition temperature of Gluio and Glnio to 
still higher values. This energy contribution also explains why the transition temperatures 
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TABLE II: Properties of the helix-coil transitions observed for the five homopolymers: Transition 
temperature Ti, half width AT, and specific heat per molecule at the transition temperature, 
C(Ti). The number in parenthesis is the error in the last digit. 



vacuum 



solvent 





Ti[K] 




C{Ti) [Kcal/mol] 


Ti[K] 


Ar"[K] 


C{Ti) [Kcal/mol] 


GluiQ 


587(14) 


161 


0.193(2) 


477(7) 


88 


0.296(3) 


GlriiQ 


584(14) 


163 


0.181(1) 


484(8) 


97 


0.268(3) 


LysiQ 


538(8) 


151 


0.216(2) 


447(10) 


98 


0.266(3) 


AsniQ 


485(19) 


193 


0.169(2) 


424(9) 


105 


0.249(4) 


AspiQ 


471(19) 


170 


0.182(2) 


415(6) 


82 


0.300(3) 




427(7) 


146 




333(2) 







Metermined at C = [C(Ti) + C{Tr. 
''The data for Ala^o are from Refs. 
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e, with T„„ from either C"(r„„) = or C"(T„ 
and are included for comparison. 



of AsriiQ and AspiQ are still higher than that of AlaiQ, although the side chains for those 
three molecules are of similar size. 

The inlay of Fig. [3] shows an interesting additional effect of side chain size. Over the 
whole temperature range the radius of gyration of Lysm is clearly larger than that of the 
other molecules. However, in the high temperature phase and in the transition regime the 
values of Tgy for Glui^ and GIuiq as well as for Asrii^ and AspiQ are very close to each 
other, despite differently sized side chains. Only in the low temperature phase do their sizes 
become recognizably different (see also below the discussion of their respective low-energy 
structures). Above this phase the size of these four polymers is apparently dominated by 
backbone fluctuations. The LysiQ side chains, however, being two or three CH2 groups longer, 
respectively, do give a considerable contribution to the size of the molecule, as expected. 

Interestingly, a second peak in the specific heat is observed for two of the homopolymers, 
GluiQ and GIuiq. The peak temperatures T2 are around 1707^, with a peak width of about 
The exact peak properties are listed in Table UTTl From the inlay of Fig. [3] it follows 
that this second peak cannot be related an additional collapse of the two molecules. We 
have argued in Ref. Il2| that this low-temperature peak is related to an ordering of the side 
chains in the form of hydrogen bonding between them. This hypothesis is supported by 
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TABLE III: Properties of the side chain ordering transitions observed for GIuiq and GIuiq : 
Transition temperature T2, half width AT, and specific heat at the transition temperature, C{T2). 
The number in parenthesis is the error in the last digit. 







vacuum 






solvent 










C{T2) [Kcal/mol] 


T2[K] 


AT''[K] 


C(T2) [Kcal/mol] 


GluiQ 


166(16) 


203 


0.176(4) 


111(10) 


166 


0.152(2) 


Glriio 


181(23) 


187 


0.161(6) 


120(17) 


146 


0.133(3) 



determined at C = [^(Ta) + C{T^^n)] /2, with T„,„ from C'(r„,„) = 0. 



Fig. m where we display the average total number n/j^ of hydrogen bonds as a function 
of temperature. Note that the fully formed helix contains six hydrogen bonds. Below 
T = 400K however, more than six hydrogen bonds are formed, with the number increasing 
with decreasing temperature. Since no hydrogen bond partners are available anymore on 
the backbone, these additional hydrogen bonds are formed between the side chains. The 
fluctuations of the hydrogen bonds, x(^) = (('^/ife ~ (^/ife))^), are shown in the inlay of Fig. HI 
In addition to the peak at Ti that appears for all polymers and that we expected from the 
helix coil transition, there is a second peak at the temperature T2 in the curves for GIuiq 
and GlriiQ. This peak corresponds to the second peak in the specific heat, and this result 
clearly indicates indeed a second transition connected to the hydrogen bond formation. It 
separates a low-temperature phase with ordered side chains from a phase at temperatures 
above T2 where only the backbone is ordered and there is only a small number of fluctuating 
side chains. The form of the side chain ordering can be seen best from the lowest energy 
conformations displayed in Fig. [61 Here, the side chains nestle along the cylinder formed by 
the helix and are stabilized by the side chain hydrogen bonds (not shown in the figure). 

Asnio and Aspio, whose side chains are only one CH2-group shorter than those of GIuiq 
and Glnio, show no such a second peak in the specific heat. This is surprising since also 
in these polymers hydrogen bonds are forming continuously with decreasing temperature 
among the side chains, like in Gluio and GIuiq (see Fig. H]). However, as the inlay of that 
figure demonstrates, this hydrogen bond formation is not accompanied by a peak in the 
hydrogen bond number fluctuations. 

We have investigated this puzzle by a more detailled inspection of the hydrogen bond for- 
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mation. Fig. [5] shows the temperature dependence of the hydrogen bonds between monomers 
that are less than four units apart in the sequence, i.e. non-backbone hydrogen bonds. It 
can be seen that in Asriio and Aspio such bonds take part in the hydrogen bond formation 
and contribute to the fluctuations already at the helix-coil transition. GIuiq and GIuiq, in 
contrast, show no such fluctuation peak. In addition, we observe that the number of these 
hydrogen bonds that are of short range in the sequence is not sufficient to account for all 
hydrogen bonds in the low temperature phase. This means that Asriio and Aspio exhibit 
some side chain hydrogen bonds in the low temperature phase that span four residues and 
are, therefore, parallel to the backbone hydrogen bonds. In contrast, practically all side 
chain hydrogen bonds in GIuiq and GIuiq cover smaller distances along the sequence. Once 
the helix is formed, residues four units apart are actually closest in space. Consequently, 
side chain hydrogen bonds in Asriio and Aspio are actually between groups closer in space 
than they are in the case of GIuiq and GIuiq. This is also suggested when comparing the 
lowest energy structures of the four polymers in Fig. [H Although in all cases the side chains 
nestle along the cylinder formed by the helix, in Asriio and Aspio they appear to be more 
stretched and densely packed. In particular, in Asriio and AspiQ the side chains appear to 
have much less possibility to find alternative hydrogen bond partners than they have in the 
case of Gluio and GlniQ, leading to this dramatic difference in the fluctuation behavior of 
both molecule groups. 

The reason for this remarkably different behavior of AsniQ and AspiQ is the apparent 
existence of a threshold length for the side chains. If these chains are not long enough, 
the number of degrees of freedom is simply too small to allow for the fluctuations observed 
for Gluio and GlriiQ. This difference can be compared to the behavior of a generic model 
for phase transitions, the Ising model^^, in one and higher dimensions. While in dimension 
d > 1 the number of effectively interacting degrees of freedom are sufficient to allow for a 
phase transition at some T > 0, they are not sufficient in the Id Ising model, leading only 
to a monotonic ordering upon decreasing the temperature to T = in that case. Similarly, 
the ordering of the side chains in AsniQ and AspiQ does not have a transition like character 
but is spread out over a wide range of temperatures. 

Such a side chain ordering behavior is also what we observe for LysiQ. Since its side 
groups could participate only in weak hydrogen bonds, they act neutral at the helix coil 
transition and finally order themselves as it is shown in Fig. El The side groups align 
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themselves along the helical axis, however less densely packed than in the case of the strong 
hydrogen bond forming side chains. We observe a minor increase of the specific heat at 
the lower temperature end of our data. However, since the error in this regime is largest, 
we consider that finding inconclusive. Rather, we have to conclude that the specific heat 
actually stabilizes over a wide range of temperatures below the helix coil transition and side 
chain ordering is actually a smooth process over this temperature range. 

We note that the structural features in the low temperature phase found here corre- 
spond to the results of theoretical structure investigations of such homopolymers in vacuum 
performed by H. Scheraga's group some 40 years ago^*^*^*^. 

So far, we have focused on the behavior of our molecules in gas phase. This research is 
in itself theoretically interesting as well as important experimentally, since the properties 
of biopolymers in gas phase have become accessible to measurements only recently21i2^. In 
particular it has been verified that helices can be stable in the gas phase up to high temper- 
atures. The high transition temperatures that we observe in our simulations are, therefore, 
not unrealistic. However, in nature proteins are usually solvated, and their function often 
depends on the details of the solvent environment. For this reason we have extended our 
investigation in a second step to that of solvated homopolymers. 

As in the case of gas phase simulations we observe for all five molecules a helix-coil 
transition characterized by peaks in the specific heat, see Fig. [3, their properties being listed 
also in Table [Tll As can be seen clearly, the width of these peaks is much narrower and 
their height is larger than in gas phase, indicating a sharper transition, a feature that has 
been found before for Polyalanine^^*^ and by us for G/uio^^. Correspondingly Fig. [8] shows 
a rapid increase in the average helicity and a sharp drop of the radius of gyration (see the 
inlay of Fig. [8]). 

All transition temperatures are shifted considerably to lower values. The reason for this 
shift is the competition between the formation of backbone hydrogen bonds that stabilize an 
a-helix, and the formation of hydrogen bonds between the backbone and the solvent in the 
coil phase, the energy contribution of the latter being described in a mean field way by the 
solvent energy term ([2]). While in vacuum the transition to the coil phase is driven solely by 
entropy, here also a part of the energy, i.e. the peptide-solvent interaction, favors the coil 
phase. These effects collaborate so that the transition takes place at a lower temperature 
and becomes sharper. However, the relative ordering of the transition temperatures among 
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the polymers remains practically the same as in the gas phase, and the explanations for this 
ordering we gave above apply here, too. 

Interestingly, again a second peak in the specific heat is observed for GIuiq and GIuiq, 
the exact properties being listed also in Table IIIII As before, this second peak cannot be 
related to a collapse of the two molecules, as can be seen from the inlay of Fig. [SI However, 
here the reason for such a behavior cannot be the formation of hydrogen bonds among the 
side chains, since they interact with the solvent instead. The lowest energy structures show 
that the side chains extend into the solvent, see Fig. [91 Also, Fig. [10] shows that the total 
number of hydrogen bonds is limited to those that stabilize the helix backbone. 

A closer inspection shows that this second peak in the specific heat is mainly due to 
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fiuctuations in the solvent contribution to the total energy. This was already found in Ref. 
for GluiQ, and we confirm it here for GlriiQ. The details of the individual contributions of 
the fiuctuations to generate the specific heat peak are rather subtle, however, and we refer 
to that reference for a more detailed discussion. We just emphasize that - despite the fact 
that the side chain ground state structures are fundamentally different from gas phase - a 
side chain ordering transition is observed in solvent for GIuiq and GIhiq. 

No such behavior can be observed for the shorter side chain polymers Asuiq and AspiQ. 
The lowest energy structures show that also in their case the side chains extend into the sol- 
vent, see Fig. [9] However, the chains are apparently just too short to generate the necessary 
amount of solvent fluctuations. Again, as before in gas phase, we encounter the situation 
that the length of these side chains appears to be just below some threshold. Surprisingly 
there is still some minor side chain hydrogen bond formation in the low temperature phase 
of AsriiQ and AspiQ, as can be seen in Fig. [TOl And as we saw in gas phase, in Fig. [TT] we 
again observe enhanced side chain hydrogen bond fluctuations at the helix-coil transition, 
albeit on a smaller level. 

Lysio, flnally, exhibits again no particular features with respect to side chain ordering. 
As to be expected, its side groups extend far into the solvent, see Fig. [9l As before in gas 
phase, this side chain ordering does not have a transition like character but is spread out 
over a wide range of temperatures. 
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IV. SUMMARY AND OUTLOOK 



We have investigated backbone and side chain ordering in five homopolymers, GIuiq, 
Glriio, Aspio, AsriiQ, and LysiQ. Those molecules were selected to cover a wide variety of 
possible interactions between the side chains of the monomers. All homopolymers undergo 
helix-coil transitions and we were able to explain the ordering of their respective transition 
temperatures. 

We found that peptides with long side chains that are capable of hydrogen bonding, 
i.e. Gluio, and GIuiq, exhibit a second transition at lower temperatures connected with side 
chain ordering. Remarkably, this occurs in gas phase as well as in solvent, i.e. independent of 
the environment, although the character of the side chain structure is different in each case. 
However, short side chains capable of hydrogen bonding, as in Aspio and AsniQ, exhibit 
no separate side chain ordering transition at temperatures below the helix-coil transition 
temperature. Instead, final side chain ordering in these polymers takes place over a wide 
temperature range, exhibiting no phase transition like character. Side chain ordering in 
Lysio, whose side groups are long and polar, also takes place over a wide temperature range 
and exhibits no phase transition like character. Again, these results are independent of the 
environment, despite the different character of the side chain structure in each case. 

Our results indicate that, the de-coupling of backbone and side-chain ordering does not 
depend on the details of the environment, but solely depend on the particular side groups 
involved. While homopolymers with some side groups show a separate ordering transition 
in both environments, some do not. Since natural proteins are heteropolymers of amino 
acids, there will be groups in a sequence that show separate side chain ordering in the 
homopolymer well as groups that do not. Whether separate ordering still occurs 

in such a situation is unclear and will depend on the details of the sequence. However, 
it is an interesting question that should be tackled in the future. Possible candidates for 
test proteins are the villin head piece subdomain HP-36 or the human parathyroid hormone 
fragment PTH(l-34) that have already been proven to be accessible by the computational 
methods used here^iiS^. 
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Figure captions: 



Fig. 1: Time series of the energy for one of the five homopolymers, GIuiq. Over the course of 
the simulation the system performs a random walk in energy. The number of indepen- 
dent visits of the low-energy region gives a measure for the quality of low temperature 
quantities. This number can be approximated by the number of tunnelings between 
preselected Ei^qk and values, -132.4 Kcal/mol and -42.6 Kcal/mol, as separately 

indicated by the two horizontal lines in the plot. 

Fig. 2: Specific heat C(T) as function of temperature T for the five homopolymers in gas 
phase as obtained from a multicanonical simulation with 5 x 10^ sweeps (10^ sweeps 
for Gluio). Error bars are included and are mostly about the symbol size or less. 

Fig. 3: Average number of helical residues, < nn > {T), as function of temperature T for 
the five homopolymers in gas phase as obtained from a multicanonical simulation with 
5 X 10^ sweeps (10^ sweeps for Gluio). The inlay shows the corresponding average 
radius of gyration < rgy > (T). Error bars are included and are mostly about the 
symbol size or less. 

Fig. 4: Average total number of hydrogen bonds, < Uhf, > (T), as function of temperature 
T for the five homopolymers in gas phase as obtained from a multicanonical simulation 
with 5 X 10^ sweeps (10^ sweeps for GIuiq). The inlay shows the corresponding fiuc- 
tuations x{T) as function of temperature T. Error bars are included and are mostly 
about the symbol size or less. 

Fig. 5: Average number of hydrogen bonds between monomers that arc less than four 
units apart in the sequence, < nf^ > (T), as function of temperature T for the five 
homopolymers in gas phase as obtained from a multicanonical simulation with 5 x 
10^ sweeps (10^ sweeps for Gluio). Note that this number is a lower limit to the 
number of non-backbone hydrogen bonds, see text. The inlay shows the corresponding 
fiuctuations x(T') as function of temperature T. Error bars are included and are mostly 
about the symbol size or less. 

Fig. 6: Lowest energy configuration of the five homopolymers in gas phase as obtained from 
a multicanonical simulation with 5 x 10^ sweeps (10^ sweeps for GIuiq) and subsequent 
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minimization. Each structure is shown in top and side view to give a better impression 
of the actual 3D structure. The pictures have been drawn with Pymol^. 

Fig. 7: Specific heat C (T) as function of temperature T for the solvated homopolymers as 
obtained from a multicanonical simulation 'with 5 x 10^ sweeps. The inlay shows the 
corresponding average radius of gyration < Tgy > (T). Error bars are included and 
are mostly about the symbol size or less. 

Fig. 8 Average number of helical residues, < uh > (T), as function of temperature T 
for the five homopolymers in gas phase as obtained from a multicanonical simulation 
with 5 X 10^ sweeps. The inlay shows the corresponding average radius of gyration 
< Tgy > (T) . Error bars are included and are mostly about the symbol size or less. 

Fig. 9: Lowest energy configuration of the five solvated homopolymers as obtained from 
a multicanonical simulation with 5 x 10^ sweeps and subsequent minimization. Each 
structure is shown in top and side view to give a better impression of the actual 3D 
structure. The pictures have been drawn with Pymol^. 

Fig. 10: Average total number of hydrogen bonds < Uhi, > (T) as function of temperature 
T for the five homopolymers in gas phase as obtained from a multicanonical simulation 
with 5 X 10^ sweeps. The inlay shows the corresponding fiuctuation x{T) as function 
of temperature T. Error bars are included and are mostly about the symbol size or 
less. 

Fig. 11: Average number of hydrogen bonds between monomers that are less than four 
units apart in the sequence, < n^^ > (T), as function of temperature T for the five 
homopolymers in gas phase as obtained from a multicanonical simulation with 5 x 10^ 
sweeps. Note that this number is a lower limit to the number of non-backbone hydrogen 
bonds, see text. The inlay shows the corresponding fiuctuation x{T) as function of 
temperature T. Error bars are included and are mostly about the symbol size or less. 

Fig. 12: TOC Graphic 
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